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FOREWARD 


Reports Hypersonic Ionizing Air Viscous Shock-Layer 
Flows Over Nonanalytic Blunt Bodies (CR-2250) and Computer 
User's Guide For a Chemically Reacting Viscous Shock-Layer 
Program (CR-2251) by Miner and Lewis should be used 
together as source or reference material. 
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INTRODUCTION 


While supersonic and hypersonic flows over blunt bodies have been of 
interest in fluid dynamics for many years, recent developments in aerodynamics 
and space flight have increasingly focused attention on the problem of predict- 
ing the blunt body flowfield. In the approach most commonly used, the flow- 
field over the body is treated in two parts, an inviscid outer flow and a vis- 
cous boundary layer. Many methods have been developed for solving the inviscid 

outer flow, as examples, the methods Inouye, Rakich and Lomax J Rizzi and 
2 3 

Inouye, and Kutler, Reinhardt and Warming. Likewise, many methods have been 

developed for solving the boundary-layer flow; two particular examples are the 

4 5 

methods of Blottner and Flugge-Lotz and Blottner. 

This approach to the problem generally worked quite well. It is, however, 
most appropriate for supersonic, high Reynolds number flows. As interest in 
hypersonic, low Reynolds number flows increased (for example, for reentry 
applications, including the space shuttle), problems were encountered in apply- 
ing first-order boundary-layer theory to such flows. Some of the problems, 
such as displacement-thickness interaction, were partially met by using second- 
order boundary-layer theory, as an example the work of Lewis . 6 Another problem 
of the boundary-layer methods is determining the edge conditions. For super- 
sonic, high Reynolds number flows, in which the boundary-layer is thin compared 
to the shock layer thickness and more specifically the entropy layer thickness, 
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it is generally adequate to consider the conditions at the boundary-layer edge 
to be the same as given by the inviscid solution at the body surface. For 
hypersonic, low Reynolds number flows in which the boundary layer is not thin, 
determining the edge conditions for the boundary layer can be most difficult 
(see, for example, Ref. 6). In the method of Blottner,^ edge conditions were 
optionally determined by tracking streamlines from the shock crossing point 
to the boundary-layer edge or by entropy-layer swallowing. 

Many of the problems (including those mentioned above) associated with 
computing viscous, hypersonic flows over blunt bodies can be overcome by the 
viscous shock-layer approach in which the entire flowfield from the body to 
the shock is treated in a unified manner. Knowledge of the shock shape is 
still required (to determine the flow properties behind the shock), but 
problems such as those of streamline tracking and displacement-thickness 
interaction are avoided. While many researchers have been involved in develop- 
ing viscous shock-layer methods, the one who achieved perhaps the greatest 

7 8 

degree of success was Davis. 

An alternative approach to obtaining solutions for hypersonic blunt body 
flows has been the use of the full Navier-Stokes equations, for example, the 
method of Jain and Adimurthy. ’ Such methods have been quite successful in 

providing solutions for the stagnation region but generally have been applied 
only about one nose radius downstream. Further, the elliptic nature of the 
equations, at least in the physical coordinates, increases the complexity of 
the solution procedure and restricts the application of the methods in the 
downstream direction. 

The first objective of the present research was to develop a method for 
predicting hypersonic, low Reynolds number flowfields over nonanalytic blunt 
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bodies with particular emphasis for the shuttle orbiter windward plane of 
symmetry. The downstream region was of considerable interest and the method 
could not be restricted to the stagnation region. A second objective was 
that the method would not be subject to the problems involved in applying 
boundary-layer theory to such flows (problems such as displacement-thickness 
interaction and streamline tracking). 

Both objectives were partially met by the viscous shock-layer methods of 
7 8 

Davis, * but his methods were restricted to analytic bodies such as hyper- 
boloids for which the pressure distribution was nearly Newtonian. Despite the 

7 8 

restriction to analytic bodies, the viscous shock-layer methods of Davis, * 
had several advantages. The principal equations were parabolic in the stream- 
wise direction, and thus there was no restriction on obtaining downstream solu- 
tions. A finite-difference method was used for solving the equations which 
gave very good accuracy in reasonably short computing times. Further, and 
quite important, for increasing Reynolds numbers the equations tend to first 
order boundary-layer equations, and thus the methods were not restricted to 
only shock-layer flow regimes but could also be applied in the boundary- layer 
regime as well. In fact, the boundary- layer equations are a subset of the 
viscous shock-layer equations. Before discussing the present work, the methods 
of Davis^ are briefly described. 

In Ref. 7, Davis developed a set of viscous shock-layer equations for a 
perfect gas valid from the body to the shock. The equations are accurate from 
the body to the shock to second order in the Reynolds number parameter, e. In 
the solution procedure used by Davis, a first global solution was obtained using 
the thin viscous shock layer (TVSL) assumption and subsequent global iterations 
were for a fully viscous shock layer (FVSL) or for TVSL. Davis considered only 
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hyperboloids and, for the first (TVSL) global iteration, used the assumption 
that the shock angle was the same as the body angle. In subsequent global 
iterations the shock angle was computed from the body angle and the previous 
global iteration value of the shock-layer thickness derivative. This tech- 
nique successfully gave the correct shock shape for the analytic bodies Davis 
considered. In Ref. 8 the governing equations were extended to treat a react- 
ing binary gas mixture. The viscous shock-layer equations were subsequently 
extended by Moss^ for nonequilibrium air (0, NO, N0 + , N, N 2 and e") and 
other gas chemistries. 

In the present work, the viscous shock-layer equations which follow the 
formulation of Davis were solved for flows over nonanalytic blunt bodies. The 
present method is for nonequilibrium, multi -component, ionizing air; disso- 
ciating oxygen is also included. Predictions of the present method were com- 
pared with the predictions of the boundary-layer method of Tong, Buckingham 
12 

and Morse for the space shuttle orbiter windward plane of symmetry at 
224,000 feet. Predictions of the present method also were compared with pre- 
dictions of perfect gas boundary-layer flow from the method given in Refs. 13 
and 14, with predictions of seven-specie, nonequilibrium boundary-layer flow 
using the method described in Refs. 15 and 16, and with the no-injection, 
experimental data of Pappas and Lee.^ Predictions of the present method were 
also compared with the results Kang and Dunn obtained with a more approxi- 
mate integral method. Predicted electron concentration profiles were compared 
with the predicted and experimental profiles given by Evans, Schexnayder, and 
Huber. 22 


4 



ANALYSIS 


In the present work, the governing equations for the viscous shock-layer 

7 8 11 

flows follow the formulation of Davis * and Moss. The shock-layer equations 

derive from the governing equations for reacting gas mixtures (such as given 

by Bird, Stewart and Lightfoot or Williams ) written for a body oriented 

coordinate system as shown in Fig. 1. The equations are first nondimension- 

alized by variables of order one at the body surface (corresponding to high 

Reynolds number, boundary-layer flows). The equations are also nondimension- 

alized by variables of order one in the outer inviscid flow (corresponding 

to the shock region). A single set of equations is then obtained by retaining 

terms from the equations in each set to second order. The resulting set of 

shock-layer equations is uniformly second-order accurate in the inverse 

Reynolds number parameter, e, from the body to the shock. Both longitudinal 

and transverse curvature are included. As given by Davis, the governing 

viscous shock-layer equations were specialized for a perfect gas^ or a binary, 

8 11 

reacting mixture of oxygen atoms and molecules. Moss gave the shock-layer 
equations for a multi -component mixture of reacting gases. 

Governing Equations 

The equations for shock-layer flows of multicomponent gases are given below. 
Continuity Equation: 

h; C( r + y cos <|)) J * pu] + |- [( 1 + <y) (r + y cos <f>) J ' P v] =0 (1) 
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s-Momentum Equation: 


1 .. 3u x au .... < L 1 

T+icy pu 3? + pV 3y + pUV 1 +Jy + ViTy 


3S 


2 a 


f au 

KU 

e 9y 

y 

( ay " 

1+Ky )_ 


2 

e y 


2k . 

j cos <t> 

au 

KU 

1+Ky 

r + y cos <p j 

, ay ' 

1+Ky 


( 2 ) 
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which becomes 
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if the thin shock-layer approximation is made. 
Energy Equation: 
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Species Conservation Equations: 
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where J.. is the diffusion mass flux term of species i, and 


Equation of State: 


MC 

poo 


With binary diffusion only and with constant binary Lewis numbers (all 
equal), the diffusion mass flux term of the species is given by 


J i “ " Pr Le i 3,y 


The species mass fractions are given by 


C i = Pj/p 

The frozen specific heat of the mixture is given by 


C n = Y C. C n 

P i—> 1 P-? 


and the mixture molecular weight is given by 


ns c -s 


The preceding equations are nondimens ional . The dimensional equations 
were nondimensionalized by the following relations: 
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The dimensionless parameters which appear in the shock-layer equations are 
given by the following relations: 


and 


_ * .. * 
Pr = C p v /k 


e = 



★ 

n 


(12a) 

(12b) 





(12c) 


For the finite-difference solution procedure, it is advantageous to 
transform the shock-layer equations. The independent and dependent variables 
(except for the species concentrations) are normalized by their local shock 
values. When the normal coordinate is normalized by the local shock-layer 
thickness, a constant number of finite-difference grid points with constant 
spacing between the body and the shock can be used. Also, interpolation is 


8 



not needed to determine the shock location and grid points in the normal 
direction need not be added. 

The transformed independent and dependent variables are 


" = y/y sh 

(13a) 

4 * S 

(13b) 

“ ' u/u sh 

(13c) 

sz 

LO 

> 

> 

II 

!> 

(13d) 

p = p/p sh 

(13e) 

P ■ P /p sh 

(13f) 

T ■ T / T sh 

(13g) 

“ * “ /lJ sh 

(13h) 

J = k/k 5h 

(131) 

r. - c „/c n 
p p p sh 

(13j) 


The transformations relating the differential expressions are 


3 _ 3 y sh 3_ 
3S " 35 ' y sh n 3n 


(14a) 


3 _ _ 1 _ 9 _ 

3y ' y sh 8 " 


and 


9 * 2 4 



(14b) 


(14c) 


where 



^sh 

d? 


(14d) 


9 


When written in the transformed 5,n coordinates, the s-momentum, energy 
and species continuity equations (Eqs. 2, 4 and 5) can be expressed in the 
following standard form for a parabolic partial differential equation: 



9W 

9n 


+ A 0 W + At + A, 


9W _ 
95 


(15) 


where W represents IF in the s momentum equation, T in the energy equation and 
C.j in the species continuity equations. The coefficients A-j through A^ are 
functions of the independent and dependent variables and may be written as 
follows: 
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Energy equation 
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Species Conservation Equation 
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In the transformed coordinates the remaining equations are 
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if the thin shock-layer approximation is made, and 


State Equation 
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The energy and species conservation equations (Eqs. 4 and 5) include the 

rate of production terms, , of species i. The w.. terms are functions of 

5 8 

both the temperature and the species concentrations. Blottner and Davis 
discuss the need for rewriting these terms so that the temperature or the 
species concentrations appear as one of the unknowns. For the energy equation, 
the production term is written so that the temperature appears as an unknown 

O 

as given by Davis as 
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where k denotes the iteration for which the solution is known and k+1 the 
iteration for which a solution is sought. Accordingly, the production term 
in the energy equation (Eq. 4) was rewritten as 


I 


i=l 


*i " *1 + T sh T * 


(23) 


and the terms w-| and ^ appear in the energy equation coefficients (Eqs. 17b 
and 17c). For the species conservation equation, the production term w?j 
written so that the species mass fractions appear as an unknown as 
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( 24 ) 


and the terms and appear in the species conservation equation coeffi- 
cients (Eqs. 18b and 18c). 

The viscous shock layer for nonequilibrium chemistry is described by 
equations (15) through (21) together with the appropriate boundary conditions 
and relations for the thermodynamic and transport properties. 


Boundary Conditions 

At the body surface, the no slip boundary conditions were imposed. For 
ri = 0, the surface conditions are 

u = 0 (25a) 

v = 0 (25b) 

and 

T - T w (25c) 

where T Is either a constant or has a specified variation. For a noncatalytic 
w 

surface, (NCW), the species boundary conditions are 



9n 


= 0 


(25d) 


The equilibrium catalytic wall (ECW) conditions are specified by 


C 1 = C. (T ) (25e) 

eq 

In the present work the surface temperatures were sufficiently low that the 
ECW condition could be approximated by a fully catalytic surface (FCW) condi- 
tion specified by 

C 0 = °* C 0 2 = °- 23456 > C N0 = ° 5 C N = °» C N0 + = 0 and C N 2 = °- 76544 (25f) 
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At the shock, the velocity components tangent and normal to the shock 
are not the same as the components tangent and normal to the body surface. 

The velocity components tangent and normal to the shock are denoted by u $h 
and v $h and the components tangent and normal to the body surface are denoted 
as u $h and v $h - The transformation relating the two sets of shock velocity 
components is 

u s h = u $h sin (a + b) + v $h cos (a + e) (26a) 

and 

v sh “ " ^sh cos + B ) + v sh sin + (26b) 

where e = tt/2 - <f>. 

For shocks of finite thickness called shock slip (SS), the shock properties 
are given by the modifed Rankine-Hugoniot relations (see Davis^’® and Cheng 2 ^) 
below. 


p sh v sh 

= - 

• sin 

a 



(27a) 

e2 "sh 1 

f au’ 

I 8 *; 

+ 

sh 

sin 

“ u sh 

= sin a cos a 

(27b) 

p sh - s1n “ "sh 

1 

Pc 

p. 

2 + 
U * 

» 00 

c-in 2 

sin a 

(27c) 

' 2 k sh 

3T| 

\dyj 

+ 

sh 

sin 

ns 

•E 

1=1 

r u sin a 

1 ‘ 2 
«> sh 

C(u sh - COS a) 2 + sin 2 a - v sh ] 



ns 





= si 

n a 
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h i 


(27d) 



i=l 

00 

oo 



and 







2 y sh 

Le i 

8C ich 



(27e) 

6 Pp sh 

9y 

— + 

sin a 

C *i - sin a C • 

sh ^ 
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The quantity p $ ^ is determined from the equation of state (6) after determining 

P sh’ T sh and C i from Eqs ‘ t 27c )> ( 27d ) and ( 27e )- In Ec l- (27d) above, the 

species enthalpy at the shock, h. , is expressed directly in terms of T 

1 sh 


sh 


before T s ^ is determined. 


With no shock slip (NSS) the Rankine-Hugoniot relations are used to 
determine the shock values. Eqs. (27a) and (27c) are unchanged. The expres- 
sions for u s ^, T g ^ and become 


U sh = COS a 


(28a) 




i=l 

and 


C i 


Ish 


- (0 


sh 


- cos 


a ) 2 /2 


+ (sin 11 


- 



■ Cl 

OO 

The shock conditions for the dependent variables (at y = 1) are 


(28b) 


(28c) 


u = v = p = P = T = l 
and 



(29a) 

(29b) 


Surface Transport 

The surface skin friction and heat transfer rates are given by the skin 
friction coefficient and Stanton number. The skin friction coefficient is given 
by 

* 

2t w 

~ ~ • * *2 (30a) 

P U 

00 OO 
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where 


T w 


* <)U 

v — *■ 
9y 


(30b) 


J w 


In terms of the nondimensionalized variables, the skin friction coefficient is 
given by 


C f = ze 


3U 

9y 


w 


(30c) 


The Stanton number is given by the expression 
* 

St = 




* * ★ ie 

p V (H - H ) 

^ 00 oo v oo W 

or in the dimensionless variables 

q„ 


St = 


n w 


H - H , 

oo W 


where 


\ = “ 


and 


0 = -e 

H W 


or 


q w = - £ 


k*<-f h^; 

8y i=l 


w 


ns 


3T 

9y 


E h i J i 


i=l 


ns 


w 


i=l 


rv T C““l 3 C • 

k ay + I fr Le i h 1 37 


(31a) 


(31b) 


(31c) 


(31d) 


(31e) 


J w 


with the restriction of constant and equal Lewis numbers. 
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THERMODYNAMIC AND TRANSPORT PROPERTIES 


The specific heat, Cp, and static enthalpy, h, are required for each of 
the species considered and for the gas mixture. Also required are the viscosity, 
y , and the thermal conductivity, k. Since the multi -component gas mixture is 
considered to be a mixture of thermally perfect gases, the thermodynamic and 
transport properties for each species were calculated using the local tempera- 
ture. The properties for the gas mixture were then determined in terms of the 
individual species properties. In this section all expressions are presented 
in terms of dimensional quantities, and the superscript star will not be used 
to denote dimensional quantities. 


Thermodynamic Properties 


The enthalpy and specific heat of the species were obtained from the 

26 29 

thermodynamic data tabulated by Browne. Browne gave tables of specific 

heat and enthalpy versus temperature in gm cal/gm mole - °K with the enthalpy 
as (H - H )/T where H was the heat of formation. In the present work the 
units were converted as 


u _ 49686 H - H* 

"i “ 1.98726 M 1 T 

and 


ft 2 /sec 2 -°R 


(32a) 


_ 49686 

p 1 1.98726 M. 


C p ; ft 2 /sec 2 -°R 


(32b) 


Second-order Lagrangian interpolation was used to obtain the values of H and 

A 

Cp from the tables. The species enthalpy and specific heat were then obtained 
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from the expressions 

h. - T H. + Ah^; ft 2 /sec 2 (33a) 

and 

C n = C n ; ft 2 /sec 2 -°R (33b) 

n n 

where Ah|j is the heat of formation of species i. The tabulated values of 
enthalpy and specific heat are given in Tables I and II and the heats of 
formation are given in Table III. 


Transport Properties 

The viscosity of each of the individual species was calculated from the 
curve fit relation 


i, - exp (C,) T k (A 1 ln T k + “l 1 ; 


qm 


cm-sec 


(34) 


29 

where , B.. and are the curve fit constants for species from Blottner 

which are given in Table IV and T^ is the local temperature in degrees Kelvin. 

2 

The units of the species viscosity were converted to lbf-sec/ft . 

The thermal conductivity of the individual species was calculated from 
the Eucken semi -empirical formula using the species viscosity and specific 
heat by the expression 


k i = M~ 


C Pi M i , 5 
R 4 


. lbf 
5 sec 6 R 


(35) 


After the viscosity and thermal conductivity of the individual species 
were calculated, the viscosity and thermal conductivity of the mixture were 
calculated using Wilke's semi -empirical relations; 
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V 


(36) 


ns 

X 

-j. 

TT 

— j. 

ns 

1=1 

r x. *. . 


j=i 

II 

X i k i 

ns 

1=1 ' 

1 E x i f u 

j=l 


Ibf-sec 



Tbf 

sec-°R 


(37) 


where = C.. M/M.. 


and 


s >. . 
ij 






-1 


In the present work, the diffusion model is limited to binary diffusion 
with the binary diffusion coefficients specified by the Lewis number from 
Eq. (12c). 


Le. = p Cp D./k 

The values of the Lewis numbers used were 1.4. 


20 


CHEMICAL REACTION MODEL 


In the present work, it is assumed that the chemical reactions proceed at 
a finite rate, and the rate of production terms, w. , of the individual species 
are needed. The production terms occur in the energy equation (Eq. 4) and the 
species conservation equations (Eq. 5). For a multi -component gas with ns 
distinct chemical species and nr simultaneous chemical reactions, the chemical 
reaction equations are written in the general stoichiometric form 


nj k nj 

L “H x i ^ i 6 ri X i 

1=1 V 1=1 


(38) 


where r = 1, 2, ... nr and nj is equal to the sum of the species and the cata- 
lytic third bodies. The quantities represent the chemical species and the 
catalytic third bodies, and the and are the stoichiometric coefficients 
for reactants and products. The rates at which the forward and backward reac- 
tions occur are specified by the forward and backward rate constants which are 
given by the equations 
C2„ 

k f = T k r exp (C0 r - Cl r /T k ) (39a) 


and 

D2„ 

\ ■ T k exp (D0 r - Dl r /T k ) 


(39b) 


where T k is the temperature in degrees Kelvin. The constants C0 r , Cl r , C2 r , 
D0 r , D1 r and D2 r depend in part on the specific reaction equations chosen. In 
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the present work reaction rate constants were matched to those used by Evans, 
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Schexnayder and Huber or by Kang and Dunn. The reaction equations and rate 
constants for these two sets of data are given in Tables V and VI. Other rate 
constants were used for test purposes (from Blottner and from Blottner ) 
and are given in Tables VII and VIII. 

With the forward and backward reaction rate constants given by Eq. (39) 
the net mass rate of production of species i per unit volume, w. , is given by 
the equation 

w nr 

r = - "VI > < L f r - L br ) < 40 > 

r=l 


where 

nj 

a r = a rj " 1 

j=l 

nj 

3 r = e nj - 1 

j=l 



a 


nj 



6 rj 


p (gm/cm^) = 0.51536 p (slugs/ft^) 


For the ns species the mass concentrations y. are given by the expressions 

J 

Cj 

Yj = jv[7 j = 1 , 2 , . . . ns 
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whereas for the catalytic third bodies the y. are given by the following 

J 

expressions 

ns 

Z (j-ns),1 Y 1 J = <ns + l), ... nj 
i=l 


The quantity Z^._ n$ j . is the third body efficiency relative to argon and is 
determined from the reaction being considered. 

As discussed previously, it is desirable to rewrite the expression for 
the rate of production terms so that the species concentrations appear as one 
of the unknowns. When rewritten in this way, the rate of production terms are 
given by the expression 


— = w? - w] C, 

p 111 


(24) 


where 


nr 


*1 = “1 Z < r H L f r + r ri L b r > 
r=l 


nr 

*1-1 + ( V Y 1 )] 

r=l 


r ri = < 


( e ri * “ a f1 -) ^ (3 rl - - a ri ) > 0 


lo 


if (e . - a .) < 0 
' p ri ri ' - 


r ri = < 


if (e ri - « r1 ) > o 


l -(3 ri - a f .) if (e ri - a r1 ) < 0 


(41a) 

(41b) 
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As discussed previously, the energy equation required the rate of produc- 
tion terms rewritten with the temperature appearing as an unknown (Eq. 22). 

That form for the rate of production term was a function of the derivative of 
Wj / p with respect to T. With temperature in degrees Kelvin, T^, the expression 
for the derivative is 


9 


p / 


M. 


nr 


= T: E <#n - «rl) [t C2 r + C V T k ' \> L f r ' < D2 r + D V T k ' ®r> V 


r=l 


(42) 


With the specification of the chemical kinetics, the system of governing 
equations for viscous shock-layer flows is complete. 

As noted above, the rate of production terms are for nonequilibrium flows. 
As conditions approach equilibrium, the present technique encounters increasing 
difficulty in obtaining a converged solution, particularly at the stagnation 

o on 

point. For dissociating oxygen, Davis following Blottner rearranged the 
rate of production terms so that equilibrium conditions could be approached 
much more closely. For dissociating oxygen, the rate of production terms 

O 

follow the procedure of Davis, and solutions may be obtained closer to equili- 
brium with the dissociating oxygen model than with the multicomponent air 
model. In fact, a lower limit in altitude or pressure currently exists below 
which the solution of the multicomponent air gas model computer code will 
not converge. The lower limit depends on the body nose radius and must be 
determined for each vehicle based upon available data. 
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METHOD OF SOLUTION 


In the present work, a finite-difference method (following Davis^) was 
used to solve the governing differential equations for the viscous shock-layer 
flows. The solutions for the continuity and n-momentum equations were obtained 
by integration with the trapezoidal rule. The s-momentum, energy and species 
conservation equations were expressed in the standard form for a parabolic 
partial differential equation 


2 

?_W. + n lH + U + A + A — =0 
a 2 A 1 9n A 2 W A 3 A 4 h 
9n 


(15) 


These equations were solved using the algorithm described by Davis. 


Solution for S-Momentum, Energy and Species Conservation Equations 
With the finite-difference grid as shown in Fig. 2, Taylor series expan- 
sions are used to relate the partial derivatives to the function values at the 
finite-difference grid points. In the £ direction the expansion for W is 

in 


W ^1 = W S1 + * 


9W 
l 35 I 


+ 0 [ ( A£ ) j 


m 


2 

Neglecting the terms of order (a?) 

W n - w n 
8W _ Vi w m 

3£ A£ 


gives the difference quotient for 8W/3£ as 

(43) 


In the n-di recti on a variable grid spacing was used and the Taylor series 
expansions are 
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w n_1 = W n - An , 

m m n- 1 


8W 1 

(4n n .-,) 2 

f 3 2 w] 

l’" 1 

+ 2! 

1 , 2 

m 

l 3n / 


+ 0 [(An n _i ) 3 ] 


m 


and 


, 1 * 1+1 ,,n , . 

W m = W m + An„ 

m m n 


' 3W 1 

n (in/ 

1 s2 “ l 

3n j 

2! 

m 

k 2 ) 


|) +0[(in n ) 3 ] 


Neglecting the terms of order (An) 3 the above equations combine to give the 
following difference quotients. 


f = a, W" +1 + b, w" + c, w "' 1 
3n I m 1 m 1 m 


(44a) 


and 


3 W - a U n+ 1 + K U n 4. r u n “l 

o ~ a o W + b, W_+ c W m 

8r| 2 2 m 2 m 2 m 

(44b) 

where 


a l = An n-l^ An n An T^ 

(45a) 

b l = K ‘ 4n n-l V(4n n 4 V 1> 

(45b) 

C 1 * ‘ 4n n /(4 Vl inj) 

(45c) 

a 2 = 2/(An n An T ) 

(45d) 

b 2 = - 2/(in„ an n _ 1 ) 

(45e) 

c 2 = 2/ (An n _i An T ) 

(45f) 
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At >T = AT1 n + An n-1 
An n ” n n+l " n n 
A Vl = \ " n n-l * 


(45g) 

(45h) 

(451) 


A more general approach is to evaluate the partial derivatives at 
(m + 0 , n)- The parameter 0 gives the following finite-difference schemes. 

/ 




0 explicit 

1/2 Crank Nicholson 

1 implicit 


The difference quotient representation of the partial derivatives in the 
n-di recti on then becomes 


aW 

9n 


■ 0 ( a i Cl + b i Ci + c i C{> + ^ 


n-1 


and 


0 ) ( a i C + b, w" + c, w"" 1 ) 

' x 1 m 1 m 1 m 

(46a) 


3 9 = 0 ( a 2 Cl + b 2 Cl + c 2 Cl) + ^ ^ 1 + b ? W " + c ? w m 1 ) 


9n 


2 "m 


Jr) n • V.O H 

2 m 2 m 


Also the function W is evaluated at (m + 0 , n) as 


w e Ci ♦ <’ - 9 > "l 


(46b) 


(46c) 


Substitution of Eqs. (43) and (46) into Eq.(15) gives the following simultaneous 
linear algebraic equations involving only W at m + 1 . 


Cl + Cl + ‘n Cl ‘ “n 
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(47) 



where n = 2, 3, . . . , N-l . The coefficients for Eq. (47) are given by 
following expressions: 

*n = < c 2 + \ c 1> 0 

B = (b« + A-. b, + A ? ) 0 + A4 /a? 
n L n 1 c x\ n 

5 n ' <*2 + V 1 ’ ® 

and 


the 


(48a) 

(48b) 

(48c) 


8 2 W 

an 2 


+ Ai 


8W 

an 


+ Ao W 


m 


(1 - 0) - A 3 + A4 


w"/A? 

m 


(48d) 


m 


m 


where A-i , A 0 , A-, and A, are the coefficients of Eq. (15) evaluated at the 
n ^n J n 4 n 

n th grid point and are given by Eqs. (16), (17) and (18) for the s-momentum, 
energy and species conservation equations. 

Assuming that 


m+1 


E W n+1 
L n w m+l 


+ F. 


(49) 


01 00 

is valid through the shock layer (see Richtmyer, also Conte and Carnahan, 
Luther and Wilkes^) then W^”j is given by 


,n-l 

Vi 


c U M 

b n-l Vl 


+ F 


n-l 


(50) 


Substituting Eq. (50) into Eq. (47) and solving for and comparing 
with Eq. (49) gives the recursion formulas 



B + A_ E 


n n-l 


(51a) 
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and 


D - A_ F •. 

F = n n n- 1 

n B n + L i , 

n n n- 1 


(51b) 


With the addition of expressions for n = 1 and n = N, the requirements 


3W 


3n 


m+1 


= 0 or W. 


1 


m+1 


W. That 
w 


for the algorithm are complete. At n = 1 , 

Eq. (47) remain valid irrespective of the finite-difference grid spacing in 
the n-di recti on requires 


F 1 = W w and E ] = 0, if W, 


1 

n+1 


W. 


w 


(52a) 


or 


F ~i ■ 0 a " d ■ i. " f 


= 0 


'm+1 


(52b) 


For n = N, the value of W is 


iN 


- W sh 


(53) 


The solution of Eq. (15) is provided by the following algorithm. Starting 
with Eqs. (52), the E n and F n are evaluated (using Eqs. (16), (17), (18), (45) 
and (48)) with n increasing from 2 to N - 1. Then the W p +1 are evaluated 
from Eq. (49) with n decreasing from N - 1 to 1 . 


Solution for Y-Momentum and Continuity Equations 

The normal momentum equation, Eq. (20a), is rewritten so that 3P/3n may 
be evaluated directly as 


w < y 


9P 


sh p sh u sh 


3n ' p sh + k y S h n) 


- p sh v sh - - 3v 


p u - 


sh 


p V 3^ 
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y sh p sh u sh v sh 
P st, <’ ' K *sh n) 


I 


(54) 


p u 


9 V , v 5v sh 

9? v sh 9? 


y sh 9V 

*sh 9? J 


with only the first term on the right side of the equation retained when the 
thin shock-layer approximation is made, Eq. (20b). With the y-momentum equation 
written in this form, Eq. (54) or (20b), the pressure derivative with respect 
to n is calculated. With F at the shock known, P sh =1 ' integration by the 
trapezoidal rule from the shock inward gives the solution of the normal momentum 
equation. 

The continuity equation, Eq. (19), when integrated yields both the normal 
velocity (v) profile and the shock-layer thickness, y ^ As given previously, 
the continuity equation is 


3_ 

95 


[y sh (r + y sh " cos p sh u sh p *3 = b; [(r + y sh " cos + > J <*sh p sh u sh p u " 


- d + k y sh n) P sh v sh p v}] 


(19) 


where j = 1 for axisymmetric flow and j = 0 for two-dimensional flow. 

The mass flux between the body (n = 0) and a given grid point n (n = n) is 
proportional to m n (with m^ denoting r\ = 1, the shock) which is given by 


.. ■ j: 


m - * *sh (r + * S h " cos p sh U sh p u dl > (55) 

Integrating Eq. (19) from 0 to n and substituting Eq. (55) gives the following 

form for the Continuity equation. 

drn 

1 d r/.. . . \.l 

‘sh p u " ■ 


dm r" , . , 

HfT J 9? [<r + * sh " cos *> ( ysh p sh u < 


11 + e *sh "> p sh v sh p v,] dn 


(56) 
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or equivalently as 


dm n . ■ 

3T = (r + y sh " cos ♦) {y sh »sh u sh 0 u " ' (1 + K y sh n) p sh v sh p v) (57) 


The term dm n /ds is obtained by evaluating Eq. (55) as s + ds/2 and s - ds/2 
and dividing by ds. The normal velocity, v, is then obtained by rearranging 
Eq. (57). 

The shock-layer thickness is obtained by integrating Eqs. (55) and (56) 
from 0 to 1 instead of from 0 to n* This gives 


m„ = 


*sh p sh u sh r p u dn + j y sh p sh u sh cos * 


r 


pun 


dn 


(58) 


and 


dn^ 

ds 


■ (r + y sh cos P)J {y sh p sh u sh - (1 y sh> p sh v sh> 


(59) 


The term dm^/ds could also be evaluated from 


% . i_ 

d5 AS 


("Vs + ds/2 " ("Vs - ds/2 


(60) 


Rearranging Eq. (60) gives 


dm. 


("Vs + ds/2 " ds + ("Vs - ds/2 


(61) 


By evaluating m^ from Eq. (61), using Eq. (59) for dm^/ds, Eq. (58) can be solved 
for the shock-layer thickness, y^. 

When written as in Eq. (19), the continuity equation is indeterminate at 
s = 0. In order to evaluate the continuity equation at the stagnation point 
the following limit expressions as c + 0 are used: 
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r -*■ 5 , cos 4> ■> 5 and u $h ->• 5 u gh where u gh = d u sh /d£. Also, at s = 0, y gh = 0. 
With these expressions Eq. 19 becomes 

+y sh ^ p sh v sh P ^ - - 0 + ’> *sh <> + *sh p sh u sh p “ (62 > 

Denoting r I a?/ 2 as r 2 anc * cos ^ A5/2 as cos ^2’ an e 9 u i va ^ ent f° rm the con " 
tinuity equation is 


-a- to + y 5h n ) J+1 Fv] - ^ 


U + 1) y sh (r 2 + y sh n cos * 2 ) J o sh u sh p u 


where p s ^ v $ ^ = - sin a = - 1 at s = 0 has been used. Integrating from 0 to 
n and rearranging terms gives the following expression for the normal velocity 
component. 

77 = I 2 1 ' (j+1) y sh p sh u sh r J r n A , „ pap ± f r 77 _ 1 


Hfej 


p (1 + y sh n) 


j+1 ' 2 


p u dn + j y sh cos 4> 2 j p u n d n J 

0 (64) 


Integrating Eq. (63) from 0 to 1 gives the following equation 

i i j+1 i 

(1+y sh )J+1= sf <j +1) *sh p sh u sh r 2 l pudri 


+ j y ch cos <t>. 


p u n dn 


which can be solved for the shock-layer thickness, y $ ^, by rearranging terms. 

An alternate method for determining the shock-layer thickness, y^, is to 
directly match the mass flux through the shock with the mass flow in the shock 
layer between the body and the shock. The mass flux through the shock, corres 


ponding to a given position on the body with radius r, is given by 


The mass flux through the shock layer is given by 

r y sh * * j 

m $ 1 = \ p u {2tt (r + y cos <|>)r dy 

-'o 


or 


* * . f . r' r 1 

">sl = U » p sh u sh hb < 2 "' J [ r ) P u <tn + j y sh cos ♦ j pun dn 


(67) 


Equating the mass flux through the shock, m s , and the mass flux through the 
shock layer, m^ , gives 


JET (r + y sh cos *) J+1 = y sh p $h u sh 


1 


1 

p~ u dn + j y sh cos <(. f pun dn 
0 J 0 

( 68 ) 


which can be rearranged to solve for the shock-layer thickness. 
Eq. (68) becomes 


1 

r + y sh cos * * y sh p sh u sh ^ M dn 

0 

and for j = 1, Eq. (68) is written as 

r 2 + 2r y sh cos * + y 2 „ cos 2 * = 2 y sh P sh u jh 

y sh cos 4 



p u dn + 


p u n dn 


For j = 0 the 
(69a) 


(69b) 


When Eq. (68) is rewritten for the limit of s = 0, the expressions for y $ ^ are 
equivalent to those obtained from Eq. (62). 
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Curvature for Nonanalytic Bodies 


In the present work, two classes of nonanalytic blunt bodies were con- 
sidered. For spherically blunted cones, the surface curvature is discontinuous 
at the sphere-cone tangent point. A continuous distribution of curvature was 
obtained by computing k from the exponential approximation to the step function. 
The approximate value of k was obtained from 

k = 1 - {1 + exp [-f (s - s tan )]} _1 (70) 


where f is a constant with a typical value of 5. Calculations were made with 

other values of f and also with k given by a true step function. These cal- 
culations showed that the effects of changing the value of f were mostly con- 
fined to the region s = s t£m + 1 and that few effects were observed for s > 4 

or 5 or for s < 0.5. 

The geometry for the second class of nonanalytic blunt bodies was speci- 
fied in tabular form. For these bodies, the curvature was calculated from 
the expression 


K 


ridM 2 + 

Id 2 z1 

1 

CM 

Jds 2 ) 

Ids 2 ) 

J 


1/2 


(71) 


or from the equivalent expression 


K = 


d 2 r 


i + 


dr 

dz 


-1.5 


(72) 


The derivatives in Eqs. (71) or (72) were evaluated with a four point walking 
least squares log-log curve fit. 
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Solution Procedure 


At each s or £ location the shock-layer equations were solved in the order 
of species, energy, s-momentum, continuity and y-momentum. At each location 
the solution was iterated until convergence was obtained for the tangential 
velocity, temperature and species concentration profiles at all points of 
the finite-difference grid. The convergence test required that 


i - w^/w* 
n n 


< 6 


where n denotes the finite-difference grid point, k denotes the previous itera- 
tion value of W n , k + 1 denotes the new iteration value of W n , W represents 
IT, T or C.. and $ is a small number, typically 0.01. After a converged solution 
was obtained at a specific location, g, the profiles were then used as initial 
profiles for obtaining a new solution at «• + A£. In this way the solution 
procedure marched downstream. 

If the governing equations were fully parabolic, only one global itera- 
tion (i.e., a solution for the entire length of the body) would be sufficient. 
However, the equations depend upon d y^/d^ (and thus the shock angle). Also, 
the y-momentum equation (in FVSL form) depends upon av/3? which is not known 
(especially at the stagnation point). The downstream dependence introduces an 
elliptic nature to the equations. The elliptic effect in the y-momentum equa- 
tion is resolved by considering TVSL flows for the first global iteration. 
Subsequent global iterations may then be FVSL using the v profiles from the 
previous global iteration. 

The elliptic effect due to d y^/ds is resolved by making a suitable 
approximation for d y^/d? for the first global iteration. Subsequent global 
iterations then use d y s ^/d? (or the corresponding shock angle) as calculated 
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from the previous global iteration. In the work of Davis * only analytic 
bodies were considered and the assumption was made for the first global itera- 
tion that the shock and body angles were the same, i.e. that d y^/d^ = 0. 

Then subsequent global iterations used the distribution of d y $ ^/d£ calculated 
from the previous global iteration. In the present work this procedure was 
also used if the body geometry was nearly analytic as was the case for the 
tabular geometry considered. However, for the spherically blunted cones, 
the pressure distribution is highly non-Newtonian and the approximation that 
the shock and body angles are equal is inappropriate. For sphere-cones, an 
initial shock shape (and thus an initial distribution of d y $ ^/dc) was deter- 
mined from a blunt body, method of characteristics procedure such as that of 
Inouye, Rakich and Lomax. ^ This shock shape was used for the first global 
iteration and for subsequent global iterations, the shock angle was calculated 
from the body angle and the smoothed distribution of d y s ^/d£ from the previous 
global iteration. 
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RESULTS AND DISCUSSION 


In the present work, the principal interest was in viscous shock-layer 

flows over nonanalytic blunt bodies such as the space shuttle orbiter and 

sphere-cones, but some predictions were made for flows over hyperboloids. 

Since for analytic bodies such as hyperboloids the present method is almost 

7 8 

identical to that of Davis, predictions of the present method for hyper- 

8 

boloids should agree almost exactly with the results of Davis. While the 
results are not presented, such was indeed the case for the 10° half-angle 

O 

hyperboloid at 225 Kft (Davis ). For this case the Reynolds number parameter 
was e = 0.197 and the Reynolds number values were Re^/ft = 8355, Re oo /R n = 690 

O 

and Re $ = 63.7. Other cases considered by Davis were from 100 Kft to 250 Kft 
with Re $ = 9555 to 23.7. The agreement of the present predictions with the 

O 

predictions of Davis for flow over a hyperboloid for this case indicates the 
accuracy of present technique. 

In the present paper, predictions are presented for viscous shock-layer 
flows over four nonanalytic blunt bodies. In the first case, the geometry 
considered was the windward plane of symmetry of a space shuttle orbiter at 
224 Kft and at a 34° angle of attack. In the second case, the body considered 
was a 20° sphere-cone at 280 and 310 Kft. This geometry was considered by 

I O pi 

Kang and Dunn as approximating the windward plane of symmetry of a space 
shuttle. The third case considered was the RAM C reentry body, a 9° sphere- 
cone with R n = 0.5 ft. at 25,000 fps at 230, 250, 265 and 275 Kft. The fourth 
case considered in the present paper was the 7.5° sphere-cone (R n = 1 in. at 
M*, = 13.4) investigated by Pappas and Lee 17 in their experiments. 
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NASA Shuttle Geometry Case 


1 2 

The case considered by Tong, Buckingham and Morse was a space shuttle 
orbiter geometry using the Rockwell Shuttle 2007 trajectory at altitudes 
between 300 and 180 Kft. They considered the windward plane of symmetry, and 
boundary-layer theory was used to predict the flow over the equivalent axi- 
symmetric body. In Ref. 12 the body geometry was specified by a series of 
polynomial curve fits which were faired into a cone of half-angle correspond- 
ing to the angle of attack. In the present work the geometry considered 
followed the approach of Tong, but difficulties were encountered in directly 
using the polynomial curve fits for the forward portion of the body. In fact, 
the surface curvature as calculated from the polynomials was not only discon- 
tinuous but also changed sign. The polynomial curve fits were replaced by a 
table of s, r and z values, and a four point walking least squares log-log 
curve fit was used to interpolate for the needed values of r and z in this 
table. This approach gave a continuous distribution of surface curvature. The 
body geometry and the corresponding shock predicted by the present method are 
shown in Fig. 3 for an entry t = 800 sec. or an altitude of 224 Kft. Predictions 
by the present method were made for this particular case, since for the space 
shuttle at 224 Kft both the first-order boundary-layer theory used by Tong and 
the viscous shock-layer theory used in the present work should be equally 
appropriate. 

The pressure distribution predicted by the present method for t = 800 sec. 
is shown in Fig. 4 with pressure distributions from Ref. 12 which were obtained 
with the tangent-cone approach. Results from the two methods agreed quite well. 

Mass fraction profiles at the stagnation point are compared in Fig. 5. 

Even though there was no attempt made to match the reaction rate data used 
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in the present method with Tong's, the present results agreed reasonably well 

with his results. The most apparent differences were in the N profiles near 

the surface. Tong's results showed a pronounced rise in the N profile from 

the wall and then a gradual decrease with a sharp decrease near the shock. The 

N profile predicted by the present method showed a slight increase and then a 

decrease similar to that of the N profile of Tong. The 0 profiles were quite 

similar. The present method predicted a slightly more rapid decrease in the 

outer portion of the flowfield. The NO profiles were also quite similar. The 

present method predicted a higher peak value of NO slightly farther from the 

shock. In considering the difference between the present results and the results 

of Tong, it should be noted that not only were the reaction rates different, but 

also the present results were FVSL and the Tong results, at the stagnation point 

only, were TVSL. Considering these differences, the agreement between the 

12 

present predictions and the predictions of Tong was quite good. 

While mass fraction profile differences did exist, there was little dif- 
ference in the predicted stagnation point heat transfer as shown in the next 
figures. Heat- transfer distributions are shown in Figs. 6 and 7 and the wall 
temperature distribution is shown in Fig. 7. The predictions of the present 
method for the noncatalytic wall (NCW) condition showed a strong dependence on 
the gas model used. Over the front of the body, the heat transfer for non- 
equilibrium air (7 species) was as much as twice that for dissociating oxygen. 
This difference decreased downstream. For the equilibrium catalytic wall (ECW) 
condition, there was little difference between the predictions of the present 
method for nonequilibrium air and dissociating oxygen. For both the ECW and 
the NCW conditions the present multi -component gas results agreed well with 
the results of Tong, but the present method did not predict the rise in heat 
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transfer at s = 1 ft. that Tong predicted when entropy-layer swallowing was 
included. It should be noted that in the results of Tong the shock shape data 
used were for a different body and that entropy-layer swallowing effects are 
strongly dependent upon the shock shape. In the present method an initial 
shock shape was assumed and then updated after each global iteration. For 
this geometry the shock shapes from the second and third global iterations 
were essentially identical. A major advantage of the viscous shock-layer 
approach over the boundary-layer approach is evident from these results. For 
the viscous shock-layer approach, the shock shape is at least partially self- 
correcting with global iteration and problems such as displacement-thickness 
interaction and shock-generated external vorticity (entropy- layer swallowing 
and boundary-layer edge conditions) do not occur. 

20° Sphere-Cone Case 

1 8-21 

Kang and Dunn considered a 20° sphere-cone with R = 4 ft. since this 
geometry reasonably approximated the windward plane of symmetry configuration 
of a space shuttle. Flowfield predictions were presented for two altitudes, 

280 and 310 Kft. Stanton number distributions predicted by the present method 
are shown in Fig. 8 with the results of Kang and Dunn. The present results 
showed a significant effect of altitude on the predicted heat-transfer distri- 
bution, but the effect was about half of that obtained by Kang and Dunn. The 
normalized heat-transfer distributions shown in Fig. 9 emphasize the differences 
between the present results and the results of Kang and Dunn. The normalized 
heat-transfer distributions showed little altitude effect for the present results 
but significant effects for the results of Kang and Dunn. Also, the shape of 
the distributions were quite different. The present results showed a sharp 
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decrease in q/q Q over the spherical cap and a pronounced change in slope near 
the sphere-cone tangent point. The results of Kang and Dunn showed a more 
gradual decrease and a gradual change in slope. 

Temperature profiles are shown in Figs. 10 and 11. The results of Kang 
and Dunn showed much less altitude dependence than did the present results. 

The present results correctly showed a thicker viscous layer at the higher 
altitude (most noticeable at s = 90) while the method of Kang and Dunn predict- 
ed almost exactly the same temperature profiles for the two altitudes, especially 
at s = 90. Moreover, the present results produced the downstream "recovery" of 
a thin boundary layer while the results of Kang and Dunn did not. 

Predictions of the present method were also made for dissociating oxygen 
(0-0 2 ) flow over the 20° sphere-cone. The present predictions for 7 species 
air agreed well with the present predictions for dissociating oxygen for heat- 
transfer distributions and temperature profiles. 

Predictions were also made (using reaction rate constants matched to those 

*1 O 0*1 

of Kang and Dunn (see Table VI) for electron concentration profiles. The 
present predictions of the electron concentration profiles are shown in Figs. 

12 and 13 with the profiles predicted by Kang and Dunn. The present predic- 
tions of N e were much lower than the predictions of Kang and Dunn. Also, the 
present method predicted a monotonic decrease of N e with increasing s; whereas, 
the method of Kang and Dunn predicted a minimum N e at s = 10 or 30 with N e at 
s = 90 higher than at s = 10, 30 and 50. Perhaps more significantly, the 
present method predicted the peak electron density in the viscous layer near 
the body; whereas, the method of Kang and Dunn predicted nearly the same pro- 
files at 310 Kft and 280 Kft while the present method predicted much fuller 
electron density profiles (especially downstream) at 310 Kft than at 280 Kft. 
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The altitude effect predicted by the present method for both temperature and 
electron density profiles seems entirely reasonable and correct. 

RAM C Reentry Case 

The RAM C flights were part of a program conducted by the NASA Langley 

Research Center for studying flow-field electron concentrations under reentry 

conditions. The body for each RAM C flight was a 9° sphere-cone with a 6 in. 

nose radius. Associated with the experimental program were theoretical studies 

18-21 

using numerical methods. For example, Kang and Dunn used a TVSL integral 
method procedure to predict electron concentration profiles for several points 
on the RAM C trajectory. Also included in the Refs. 18-21 were other flow-field 
quantities such as temperature profiles and surface heat-transfer distributions. 

The results presented by Kang and Dunn were for the higher altitude points 
on the RAM C trajectory where viscous shock-layer theory would be more appro- 
priate. Predictions of the RAM C electron concentrations have also been made 

22 34 

by other researchers. For example, Evans, Schexnayder and Huber * applied 
two different boundary-layer methods and obtained reasonable agreement with the 
experimental data with both methods. The use of the boundary-layer approach 
limited Evans et al . to consideration of the lower altitude points of the RAM C 
trajectory. However, the 230 Kft point was considered by Evans et al . and by 
Kang and Dunn. For the 230,000 ft. altitude point, predictions of the present 
finite-difference, viscous shock-layer method were compared with the predic- 
tions of a finite-difference, boundary-layer method (Evans, Schexnayder and 
Huber ) and the predictions of the TVSL integral method of Kang and Dunn. 

In the present work the principal emphasis was not on predicting electron con- 
centrations but rather was on predicting the hypersonic, viscous flowfield 
over nonanalytic blunt bodies with electron concentrations a part of the 
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flowfield predictions. The present results for the RAM C at 230 Kft included 
viscous shock-layer predictions for three gas models (perfect gas, dissociating 
oxygen, and multi -component, ionizing air). Also included were results of the 
inviscid, perfect gas method of Inouye, Rakich and LomaxJ a perfect gas, 
boundary-layer method (Refs. 13 and 14) and a seven species, nonequilibrium, 
ionizing air boundary-layer method (Refs. 14 and 15). Some predictions for 
flow over a 9° half-angle hyperboloid were included in addition to the predic- 
tions for the 9° sphere cone. 

One reason for emphasizing the RAM C conditions at 230 Kft was experimental 

data and other numerical results for distributions of Stanton number, temperature 

behind the shock and shock-layer thickness in addition to temperature and 

electron concentration profiles. The experimental data and the numerical re- 

22 

suits of Evans, Schexnayder and Huber were for electron concentration profiles 

only. The differences between the present predictions and the numerical results 
1 8-21 

of Kang and Dunn were unexpectedly large. Most of the discussion of the 
differences between the present results and the results of Kang and Dunn will 
be deferred to the end of this section. 

Distributions of temperature behind the shock are shown in Fig. 14 for 
the RAM C sphere-cone and a 9° half-angle hyperboloid. Although the end of 
RAM C body was at s = 9.2 the results given by Kang and Dunn went to s = 90 and 
the present predictions were extended to s = 100. At the stagnation point there 
were small differences in T s ^ due to differences in the gas model but the sphere- 
cone and hyperboloid gave the same value for the same gas model. The differences 
in T $ ^ due to differences in gas model for the present results were greatest 
at the stagnation point and decreased as the shock became more oblique. In 
the downstream portion of the flow, the principal differences in T s ^ were due 
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to the differences in the bodies. The present results showed the expected 
differences in T ^ distributions for hyperboloids and sphere-cones and correctly 
predicted the distributions of T ^ coming together at s - 80. In contrast to 
the present results, Kang and Dunn obtained a quite different distribution 

° f V 

Distributions of shock-layer thickness are compared in Fig. 15. While the 
present results show some distinct differences, they were in substantial agree- 
ment, especially for s > 20. The principal differences in the present results 
were again due to the differences in the bodies considered. The present shock- 
layer predictions showed only small differences in y $ ^ due to the chemistry. 

For the sphere-cone, the viscous shock-layer results were in reasonable agree- 
ment with the inviscid results but did not show as sharp a decrease in y $ ^ 
between s = 15 and 40. The present results for the hyperboloid were lower than 
the sphere-cone results for s < 20 but were about 50% greater for s > 50. 

Temperature profiles are shown in Figs. 16 and 17. The profiles in Fig. 

16 were for the probe location. The present method profile was for s = 8.8, 
the probe location. The profile of Kang and Dunn (for s = 10) was the profile 
closest to the probe location which they gave (the next closest profiles were 
for s = 3,0 and 20.0). The profiles in Fig. 17 were for s = 90. The profiles 
for the hyperboloid showed a smoother transition from the outer flow to an 
inner, viscous flow than did the present sphere-cone profiles. The profiles 
for dissociating oxygen showed a ten percent lower peak temperature than did 
the profiles for ionizing air. Despite these small differences, the four 
present profiles showed a very distinct inner viscous region (y/y s ^ = 0.0 to 
0.5) and an outer inviscid region (y/y sh = 0.5 to 1.0). While not shown, the 
velocity profiles also indicated the edge of the viscous layer at y/y s ^ ' 0.5. 
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The present viscous shock-layer temperature profiles have also been compared 

(not shown here) with the temperature profiles predicted by the seven species 

15 16 

boundary-layer method * for flow over the 9° sphere-cone. The boundary-layer 
and viscous shock-layer profiles showed differences in peak and edge temperatures 
but were quite similar in character and the boundary-layer profiles also indicated 
the edge of the viscous layer at y - 1.0 as did the viscous shock-layer profiles. 
Further, a comparison of the present profiles, as shown in Figs. 16 and 17, 
clearly showed the distinct downstream development of an outer inviscid flow 
and an inner viscous flow. 

Electron concentration profiles for the RAM C at 230 Kft are shown in 

Figs. 18 and 19. The data that the present results are compared with were taken 

from figures in Ref. 22. In Fig. 18, the present results are compared with the 

22 

experimental data and the results Evans, Schexnayder and Huber obtained with 

30 

a very reliable finite-difference, boundary-! aye r method (i.e. Blottner ). 

The reaction rate constants used for the present results in Fig. 18 were matched 
to those of Ref. 22 (Table V). The present results agreed reasonably well with 
the experimental data and with the boundary-layer theory results of Ref. 22 
as to level of ionization and quite well with the boundary-layer theory results 
as to character of the N e profile. The present viscous shock-layer theory 
predicted a higher temperature in the viscous layer than did the boundary-layer 
theory of Ref. 22, and this difference accounted for most of the difference in 
the N e profiles. 

Predictions of electron concentrations were also influenced by the reaction 
rates used. Predictions were made for the RAM C case using reaction rate 
constants matched to those of Kang and Dunn^ (Table VI). The profiles pre- 
dicted using the two different sets of rate constants are compared in Fig. 19. 


45 


The principal difference was a one-third reduction in peak N g using Ref. 19 

19 

rate constants. The experimental data and the TVSL results of Kang and Dunn 

were taken from Ref. 22. Except for the experimental data between y = 10 

and 14 cm, which were affected by probe heating, the use of the Ref. 19 

reaction rate constants improved the agreement between the predictions of the 

present method and the experimental data, at least for y < 10 cm. However, as 

22 

noted by Evans, Schexnayder and Huber, the experimental data did not support 

1 9 

the upswing in N e near the shock that the Kang and Dunn results gave. The 
present results also did not show such an upswing, but rather showed the 
opposite trend. 

Heat-transfer distributions are shown in Fig. 20. The present results 
include boundary-layer (perfect gas^ 3 ’^ 4 and seven species nonequilibrium 
air^’16) and viscous shock-layer (perfect gas, dissociating oxygen, and seven 
species nonequilibrium air) predictions for the RAM C 9° sphere-cone and seven 
specie TVSL predictions for a 9° half-angle hyperboloid. All of the present 
results were in good agreement for this case. Some differences due to chemistry 
and geometry did exist, but the agreement was good despite the diversity of 
methods, chemistry and geometry. 

The heat-transfer rates at the stagnation point are not clearly shown in 
Fig. 20 but are given in Table IX. The present results for TVSL, seven species, 
1st iteration were obtained using the method of Ref. 16. The present results 
showed some distinct differences but agreed reasonably well, at least compared 
with the results of Kang. 

In Refs. 18-21, Kang and Dunn presented temperature profiles, heat transfer 
rates, shock temperature and shock-layer thickness distributions for the RAM C 
reentry body only at 230 Kft. However, predictions of electron concentration 
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profiles for s = 8.8 were given for 250, 265 and 275 Kft as well as for 230 Kft. 
In Ref. 21, Kang and Dunn noted that at 230 Kft the electron concentration 
profiles predicted with a single ionizing species (seven total species) and 
with five ionizing species (eleven total species) were essentially the same, 
but that at 275 Kft the five ionizing species chemistry model predicted electron 
concentrations an order of magnitude higher than the single ionizing species 
(N0 + ) model. For these two altitudes they gave N e profiles predicted with both 
7 and 11 species, but at 250 and 265 Kft they presented N e profiles predicted 
with the 11 species model only. 

In Figs. 21 and 22, electron concentration profiles predicted by the 

present method are compared with the experimental data (from Ref. 22) and the 

18-21 

results of Kang and Dunn " . The present results are for a single ionizing 

species (7 species total). In Fig. 21 the present results are for no shock 
slip (NSS) and the results of Kang and Dunn are for 7 species. In Fig. 22, 
the present results are for shock slip (SS) and the results of Kang and Dunn 
are for 11 species. 

Comparing the present predictions without shock slip with Kang and Dunn's 

predictions with shock slip, Fig. 21, the present predictions agreed reasonably 

18-21 

well as to level of ionization with the predictions of Kang and Dunn 
with the N0 + only model. The agreement between the present predictions and 
the experimental data was reasonably good at 230, 250 and 265 Kft, but at 
275 Kft the present method, without shock slip, significantly underpredicted 
N e . In contrast, the present predictions of electron concentration profiles, 
when shock slip was included , as shown in Fig. 22, agreed quite well with the 
experimental data, even at 275 Kft. Without shock slip, the species concen- 
trations behind the shock were the same as in the free stream and thus 
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C N0 + = c e - = ^Ith s * 10c * ( S ^ 1 P> however, a finite concentration of N0 + and 
thus e" was permitted behind the shock and diffusion carried the ions to the 
shock zone. While the electron density at the shock was quite low with shock 
slip, the nonzero value of 0^+ behind the shock raised the electron density 
profile, especially at the higher altitudes. Also at 275 Kft, with shock slip 
a thicker viscous shock layer was predicted, as was an increased static tempera- 
ture and a decreased density distribution from the maximum values in the layer 
to the shock. Thus the present method predicted an increased ionization due to 
the higher temperature and decreased deionization due to the lower density. 

Further, the present predictions showed correctly two trends not shown by 
the predictions of Kang and Dunn: (i) the viscous-layer thickness increased 

with altitude, and (ii) the present predictions showed that the peak of the 
electron concentration profile occurred within the viscous layer and not at the 
shock. In Ref. 21, the explanation given by Kang and Dunn for the peak of the 
N e profile occurring at the shock was that the temperature immediately behind 
the shock was quite low, with an accompanying high density. Thus, it was 
reasonable for the peak electron concentration to occur at the shock since, 
with shock slip, there was a finite electron concentration at the shock. If 
this were modified to state that the peak occurred near tne shock, this would 
appear to be reasonable. In fact, such results were indicated by the present 
predictions at 265 and 275 Kft. However, it does not appear reasonable that 
the electron concentration profiles would show a peak at the shock, or such a 

TO p“| 

strong upswing toward the shock as predicted by Kang and Dunn. 

In the preceding discussion, the differences between the present results 

18-21 

and the results of Kang and Dunn were only briefly mentioned. In Refs. 
18-21, a large number of results were presented but comparisons were made with 
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no other numerical results, and the only experimental results with which 

comparisons were made were for electron concentration profiles. As mentioned 

above, their results apparently agreed well with the N e profiles but as Evans, 

22 

Schexnayder and Huber noted the experimental data did not support the upswing 
in the N e profile near the shock as obtained by Kang and Dunn. It was also 
observed in Ref. 22 that the method of Kang and Dunn overpredicted the N e 
values measured by the microwave reflectometers on the RAM C-II at the more 
forward body stations (s =0.8 and 2.1) by factors as large as 20; whereas for 
233 Kft and lower, the results of the boundary-layer theory used in Ref. 22 
agreed with the reflectometer data at all body stations. These two notes from 
Ref. 22 raise questions about the results given in Refs. 18-21. Further 
questions must be raised by the differences between the results from Refs. 18-21 
and the present results. 

As mentioned above the emphasis of the present work was on predicting the 

hypersonic viscous flowfield over nonanalytic blunt bodies including spherically 

blunted cones with electron concentration profiles only a part of the flowfield 

predictions. Since electron concentration profiles are subject to changes in 

reaction rate constants as well as changes in temperature profiles, mean flow- 

field quantities such as heat-transfer distributions would be a more reliable 

measure of method accuracy. For example, predictions for the RAM C at 230 Kft 

2Q 

were also made with the reaction equations and rate constants from Blottner 
(Table VII). Predictions of heat-transfer rates were only slightly affected, 
but the predictions of the electron concentrations were an order of magnitude 
lower than when other rate data (Tables V, VI and VIII) were used. Since there 
were no experimental data -For y s ^ or St distributions, calculations were made 
using well known, well established, independent methods (inviscid flowfield 
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technique of Ref. 1 and boundary-layer flows from Refs. 13-16). The predic- 
tions of the present viscous shock-layer methods agreed well with the results 

of these independent methods. In contrast, there were large differences 

1 8-21 

between the present predictions and the results of Kang and Dunn. 

The value of T $h at s = 0 obtained by Kang and Dunn as shown in Fig. 14 
was nearly the same as for the present results, but for s > 1 their values of 
T gh were distinctly lower than the present T^ values. For s > 80 their 
value of T $h was only about 60% of the present values. The distributions of 
T s ^ obtained by Kang and Dunn should imply a lower shock angle (for s > 5 or 
10) than that in the present results. However, the shock-layer thickness dis- 
tributions shown in Fig. 15 seem to clearly indicate that (for s > 10) the 
shock angle obtained by Kang and Dunn was considerably greater than that in 
the present results. In fact, at s = 10 the values of y $ ^ were all approximately 
unity, but at s = 90 the present results gave y $ ^ - 1.9 for the inviscid flow 
over the sphere-cone, y $ ^ = 2.2 for the sphere-cone shock-layer flows and 
*sh 1 3 - 5 for the hyperboloid shock-layer flows; whereas Kang and Dunn's results 
gave y s ^ - 9. Thus, the trends of the T $ ^ and y $h distributions from Refs. 

18-21 were distinctly contradictory and inconsistent. 

The temperature profiles as shown in Figs. 16 and 17 also differed markedly. 
The profiles of Kang and Dunn were nearly the same shape at s = 10 as at s = 90, 
while the present profiles showed a strong downstream influence. The present 
profiles showed a distinct outer inviscid region, somewhat weak at s = 8.8 but 
quite clear at s = 90. The present profiles showed about the same peak value 
of T at s = 8.8 and 90. In the profiles from Refs. 18-21, the peak value 
increased from T = 7200° K at s = 10 to T - 8200° K at s = 90. 
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With the temperature profiles as different as shown in Fig. 16, differences 
in electron concentration profiles were expected, but not the differences shown 
by the profiles in Fig. 19. The present results gave a peak value of N e at 
y = 5 cm while the peak of the temperature profile was at y - 3.5 cm, but T at 

y = 5 cm was only slightly lower than T at y = 3.5 cm. This was quite reasonable, 

but in the results of Kang and Dunn, the peak value of N e was at the shock 

(y - 1 7 cm) while the peak value of T was at y = 9 cm. Further, T at the shock, 

where the peak in N e occurred, was less than one-eighth of the peak temperature. 
Also, while the N e profile of Kang and Dunn apparently agreed well with the 
experimental data, the peak value of N e was two or three times the peak of the 
experimental data (points affected by probe heating excluded), twice the peak 

of the present results and three times the peak obtained by Evans, Schexnayder 

22 

and Huber. 

At the stagnation point, the heat transfer obtained by Kang and Dunn was 
one-half to one-third of that predicted by the finite-difference methods used 
in the present work (see Table IX). In contrast to the stagnation point results, 
over most of the conical portion of the body the heat transfer predicted by Kang 
and Dunn was two to six times that predicted by the boundary-layer and viscous 
shock-layer methods used in the present work. However, the temperature profiles 
shown in Figs. 16 and 17 would apparently indicate that the heat transfer 
obtained by Kang and Dunn should have been lower than the present results. 

Ames Experimental Case 

A measure of the validity of a theory is the agreement with experimental 
data. For the shuttle configuration, flight heat-transfer data are some years 
in the future and, in general, wind-tunnel data for shuttle configurations 
are not readily available outside of the NASA and some contractors. One set 
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of experimental hypersonic wind tunnel data which has been published is that of 

Pappas and Lee^ at the NASA Ames Research Center for flow over a 7.5° sphere- 

cone with R n = 1 in. In the experimental program, surface pressure and heat- 

transfer distributions were measured at Mach 13 with varying rates of injection 

of foreign gases. Included in the experimental data were distributions for the 

no injection case. Experimental and present predictions of pressure and heat- 

transfer distributions are shown in Fig. 23 and 24. Also shown in these figures 

1 5 

is the previous first-order boundary-layer theory of Lewis, Adams and Gilley 
including transverse curvature and displacement-thickness interaction for the 
Ames conditions. The results from Ref. 15 were obtained using a global itera- 
tion for determining the displacement-thickness interaction effects, and the 
inviscid body pressure for the effective body was obtained using a blunt body, 
method of characteristics procedure similar to that of Ref. 1. The present 
theory did not compare as well with the experimental data as did the previous 
boundary-layer with viscous interaction included. In the present viscous 
shock-layer method, the effect of the discontinuity in surface curvature, k, 
was most distinct immediately upstream of the sphere-cone tangent point and 
for a short distance downstream. The sphere-cone considered by Pappas and 
Lee^ ended at s - 5, and almost all of this body was within the length affect- 
ed by the discontinuity in k. Despite the effect of the discontinuity in <, 
the agreement between the experimental data and the predictions of the present 
viscous shock-layer theory was quite good. 

While the RAM C, 230 Kft, conditions were quite different from the Ames 

conditions, the Reynolds numbers were of the same order (Re oo /R n = 4315 for the 

RAM C conditions and Re ro /R n = 1515 for the Ames conditions). The shock Reynolds 

numbers were also similar (RAM C, Re = 269; Ames, Re = 193) and the values 

s s 


52 



of the Reynolds number parameter were nearly the same (RAM C, e = 0.0965; 

Ames, e = 0.0980). The Reynolds number similarity between the two cases 

should allow comparison of the normalized heat-transfer distributions. Fig. 

25 shows the same data as Fig. 24 but with the present results for the RAM C 

conditions added. The present predictions for the RAM C and Ames conditions 

were in quite good agreement even though there was a difference in cone angle 

(and thus in the location of the sphere-cone tangent points). Further, the 

present viscous shock-layer results for the RAM C conditions agreed well with 

18 21 

the Ames experimental data. The results of Kang and Dunn for the RAM C 
are also shown in Fig. 25 in normalized form. Fig. 25 clearly shows that for 
s > 3 the results of Kang and Dunn were higher by an order of magnitude or more 
than the present results. A comparison of the results of Kang and Dunn for the 
RAM C with the Ames experimental data showed a difference by a factor of 11 or 
12 at s = 4 or 4.5. The values of Re^/Rp and Re $ given above indicate that 
the Ames conditions were at least as much in a viscous shock-layer regime as 
the RAM C, 230 Kft, conditions and it is most surprising that the trend of the 
results of Kang and Dunn did not agree better with the experimental data 
of Pappas and Lee.^ 


Computing Time Required 

Some of the computing times required for the RAM C conditions are given 
in Table X. These computing times were obtained on the IBM 370/158 system of 
the Computing Center of the Virginia Polytechnic Institute and State University. 
The inviscid gas model calculation used the blunt body, method of characteristics 
technique of Ref. 1. The PG boundary-layer (BL) calculation was made with the 
method of Refs. 13, 14. The viscous shock-layer (VSL) computing times were 
for the present method. The inviscid method generated the pressure distribution 
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for the BL calculation and the initial shock shape for the VSL calculations. 

In the VSL method, the first global iteration was for O-O 2 and subsequent 
global iterations were for either O-O 2 or the 7 sp gas model. 

The computing time (not shown) for a PG VSL was nearly the same as for 
the PG BL. As shown, a global iteration for the O-O 2 VSL required about twice 
the computing time of the PG BL. For O-Og* when the second global iteration 
was TVSL, the computing time was almost the same as for the first global itera- 
tion. However, when the second iteration was FVSL, the second iteration 
required three times the computing time of the first iteration when the same 
value of N was used, and twice the computing time when the step size restriction 
was relaxed. 

For the 7 sp gas model, the step size restriction was relaxed. The TVSL 
global iteration required six times the computing time required for O-Og, and 
for the FVSL global iteration, the computing time was 3.5 to 5 times that 
required for the O-O 2 global iteration. 

These computing times show that if the differences between the results 
for TVSL and FVSL are not great, a significant amount of computing time can 
be saved by using the TVSL model. Further, unless the 7 sp model is required, 
the computing time can be greatly reduced by using the O-O 2 gas model. 
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CONCLUSIONS 


The results of the present finite-difference method for predicting hyper- 
sonic viscous shock-layer flows over nonanalytic blunt bodies were compared 
with predictions of other finite-difference methods and with experimental data. 
For the windward plane of symmetry of a shuttle orbiter configuration, the 
predictions of the present method agreed well with the boundary-layer predic- 
tions of Tong, Buckingham and Morse. For the low Reynolds number flows over 
a 20° sphere-cone, a pseudo-shuttle configuration, the predictions of the 
present method appeared quite reasonable. Also, the altitude effects on the 
temperature and electron concentration profiles were correctly predicted by 
the present method as was the "recovery" of a thin boundary layer on the down- 
stream portion of the sphere-cone. 

The agreement of the prediction of the present method with experimental 
data further tends to verify the appropriateness of the present method. For 
the RAM C, the predictions of the present method agreed well with the experi- 
mental electron concentration profiles, with the boundary-layer predictions 
of Evans, Schexnayder and Huber and with other boundary-layer predictions of 
heat-transfer rate distributions (for perfect gas and nonequilibrium air). 
Predictions of the present method also agreed well with the experimental data 
of Pappas and Lee for pressure and heat-transfer rate distributions. 

The present viscous shock-layer method, accurate to second order in the 
Reynolds number parameter e, eliminates most of the problems encountered in 
applying boundary-layer theory to hypersonic, low Reynolds number flows over 
nonanalytic blunt bodies. The comparisons of predictions of the present 
method with the results of Kang and Dunn indicate that the present method 
is clearly superior to the more approximate method of Kang and Dunn. 
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TABLE I. Species Enthalpy vs. Temperature 
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TABLE II. Species Specific Heat vs. Temperature 
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TABLE III. Species Heats of Formation 
and Molecular Weights 


Species 

a h*j (ft 2 /sec 2 ) 

H* (Kcal/mole) 


0 

1.661 x 

10 8 

58.9725 

16.000 

°2 

0 


0 

32.000 

NO 

3.225 x 

10 7 

21 .477 

30.008 

N 

3.619 x 

10 8 

112.507 

14.008 

N0 + 

3.5341 


235.836 

30.008 

n 2 

0 


0 

28.016 



10.388 x 10 8 u * 
23.053 M. n 
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TABLE IV. Species Viscosity Curve Fit Constants 


* C. T (A. In T. + B.) , . , 

y . = exp l T k v i k 1 ' (gm/cm-sec) 


T^, degrees Kelvin 


Species 

A i 

B i 

c i 

0 

0.019558 

0.438511 

-11 .6235 

°2 

0.038271 

0.021076 

- 9.5989 

NO 

0.042501 

-0.018874 

- 9.6197 

N 

0.0085863 

0.6463 

-12.581 

N0 + 

0.042501 

-0.018874 

- 9.6197 

Np 

0.048349 

-0.022485 

- 9.9827 
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TABLE V. Reaction Equations and Reaction 
Rate Constants from Evans, Schexnayder 
and Huber (Ref. 22) 


Reaction Equations 


r = 1 

0 2 + M 1 

■> 

4- 

i 

s" 

1 

4 - 

1 

O 

1 CSJ 

1 


2 

°2 t0 2 

-> 

4- 

20 + 0 2 


3 

°2 + N 2 

-* 

4- 

20 + N 2 


4 

o 

+ 

(M 

O 

-* 

4- 

20 + 0 


5 

N 2 + M 2 

-* 

4- 

2N + M 2 


6 

N 2 + N 

-> 

•4- 

2N + N 


7 

N 2 + N 2 

t 

2N + N 2 


8 

NO + M 3 

4- 

N + 0 

+ M 3 

9 

NO + M 4 

-V 

4- 

N + 0 

+ m 4 

10 

NO + 0 

->- 

4- 

0 2 + N 





NO + N 


11 

n 2 + 0 

4- 


12 

N + 0 

-> 

4- 

N0 + + e" 



Matrix of Catalytic Third Bodies; Z^._ ns ^ ^ 


0 0 2 NO N N0 + N 2 




i=1 

2 

3 

4 

5 

6 

(j-ns) = 1 

M 1 

0 

0 

1 

1 

0 

0 

2 

m 2 

1 

1 

1 

0 

0 

0 

3 

M 3 

0 

1 

0 

0 

0 

1 

4 

M 4 

1 

0 

1 

1 

0 

0 

5 

e" 

0 

0 

0 

0 

1 

0 
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TABLE VI. Reaction Equations and Reaction Rate 
Constants from Kang and Dunn (Ref. 21) 


Reaction Equations 


r = 1 

°2 + M 1 


20 

+ M-| 


2 

CM 

O 

+ 

CM 

O 

-► 

4- 

20 

+ °2 


3 

°2 + N 2 

-V 

4 

20 

+ N 2 


4 

°2 + 0 

-* 

4 - 

20 

+ 0 


5 

N 2 + M 2 

4 

4 

2N 

+ M 2 


6 

N 2 + N 

-> 

4 

2N 

+ N 


7 

n 2 + n 2 

4 

4 

2N 

+ N 2 


8 

NO + M 3 

4r 

N 

+ 0 

+ m 3 

9 

NO + M 4 

4 - 

N 

+ 0 

+ m 4 

10 

NO + 0 


CM 

O 

+ N 


11 

n 2 + 0 

-4 

4- 

NO 

+ N 


12 

N + 0 


N0 + 

+ e“ 



Matrix of Catalytic Third Bodies; 




0 

°2 

NO 

N 

N0 + 

n 2 



i=l 

2 

3 

4 

5 

6 

(j-ns) = 1 

M 1 

0 

0 

1 

1 

0 

0 

2 

m 2 

1 

1 

1 

0 

0 

0 

3 

M 3 

0 

1 

0 

0 

0 

1 

4 

M 4 

1 

0 

1 

1 

0 

0 

5 

e" 

0 

0 

0 

0 

1 

0 
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TABLE VII. Reaction Equations and Reaction 
Rate Constants from Blottner (Ref. 29) 


Reaction Equations 


1 

1 

C\1 

1 O 

1 + 

1 CM 

O 

1 

*4* 

20 

i 

CM 

1 O 

1 + 

2 

o 

+ 

CM 

O 

-4- 

20 

+ 0 

3 

°2 + M l 

-4* 

-4- 

20 

+ M 1 

4 

n 2 + n 2 

-4 

-4- 

2N 

+ N 2 

5 

N 2 + N 

-4 

- 4 - 

2N 

+ N 

6 

n 2 + m 2 

-4 

-4- 

2N 

+ m 2 

7 

NO + M 3 

-4 

- 4 - 

N 

+ 0 

8 

NO + 0 

-4 

-4- 

°2 

+ N 

9 

n 2 + o 

-4- 

NO 

+ N 

10 

n 2 + o 2 

-4 

-4- 


2 NO 

11 

N + 0 

-4 

• 4 “ 

N0 + 

+ e - 


Matrix of Catalytic Third Bodies; Z^_ ns ^ ^ 




0 

i-1 

°2 

2 

NO 

3 

N 

4 

N0 + 

5 

N 2 

6 

(j-ns) = 1 

M 1 

0 

0 

1 

1 

0 

1 

2 

m 2 

1 

1 

1 

0 

0 

0 

3 

M 3 

1 

1 

0 

1 

0 

1 

4 

M 4 

0 

0 

0 

0 

1 

0 
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TABLE VIII. Reaction Equations and Reaction 
Rate Constants from Blottner (Ref. 30) 


Reaction Equations 


r = 1 

°2 + M 1 

-*■ 

20 

+ M 1 

2 

n 2 + m 2 

-y 

+■ 

2N 

+ m 2 

3 

N 2 + N 

-> 

2N 

+ N 

4 

NO + M 3 

-y 

4 - 

N 

+ 0 

5 

NO + 0 

-y 

4 - 

°2 

+ N 

6 

n 2 + 0 

-v 

«- 

NO 

+ N 

7 

N + 0 

4- 

N0 + 

+ e' 


Catalytic Third Bodies Efficiencies Relative 




to Argon; Z 

(j-ns) 






O 

O 

ro 

NO 

N 

N0 + 

N 2 



i-i 

2 

3 

4 

5 

6 

(j-ns) = 1 

M 1 

25 

9 

1 

1 

0 

2 

2 

m 2 

1 

1 

1 

0 

0 

2.5 

3 

M 3 

20 

1 

20 

20 

0 

1 

4 

e’ 

0 

0 

0 

0 

1 

0 
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22.9238 9.03 x 10 32,400 1/2 44.3369 1.80 


TABLE IX. Stagnation Heat Transfer for 


9° Sphere-Cone, R n = 6 in., RAM C Conditions, 230 Kft 


Model 

Gas 

Iter 

Wall 

Shock 

q, BTU/ft 2 -sec 

BL 

PG 


Present Results 


231 .074 

BL 

7 sp 

— 

ECW 

— 

177.250 

FVSL 

PG 

3 

— 

SS 

123.973 

FVSL 

CVJ 

o 

1 

o 

3 

ECW 

SS 

238.742 

TVSL 

CM 

O 

1 

o 

2 

ECW 

NSS 

190.278 

TVSL 

7 sp 

2 

ECW 

NSS 

156.883 

TVSL 

7 sp 

1 

ECW 

SS 

252.539 

TVSL 

11 sp 

— 

Kang Results 
ECW 

SS 

87.772 
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Table X. Computing Times for RAM C 
Sphere-Cone to s/R n = 120 a 


Global Computing 

Iteration Gas Viscous No. of No. of Time e 

Number** Model Model Stations Iteract1ons c N d Min: sec 


- 

PG 

Invlscld 

— 

— 

- 

5:35 

- 

PG 

BL 

61 

189 

3 

1:34 

1st 

CVI 

O 

1 

o 

TVSL 

52 

271 

3 

3:29 

2nd 

o-o 2 

TVSL 

58 

291 

3 

3:34 

2nd 

0-0 2 

FVSL 

100 

715 

3 

10:00 

2nd 

CM 

o 

1 

o 

FVSL 

62 

474 

4 

6:58 

2nd 

7 sp 

TVSL 

42 

241 

4 

20:40 

2nd 

7 sp 

FVSL 

57 

415 

4 

35:45 


^Convergence test of 1% for velocity, temperature and species profiles at each 
grid point. 

^Data are for the Indicated global Iteration only. 

c Total number of station Iterations for the global Iteration. 

d The s step size was doubled If a converged solution was obtained with the 
number of station Iterations <_ N. 

Execution time; IBM 370/158. 
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Figure 1. Coordinate System for Viscous Shock-Layer Over Blunt Bodies 
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Figu 




nO 

CT\ 



Figure 3. Body and Shock Geometry for NASA Shuttle-like Body 
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Mass Fraction 



Figure 5. Stagnation Point Species Mass Fractions 
for NASA Shuttle-like Body 
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Figure 6. Heat-Transfer Distributions for NASA Shuttle-like Body 
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* » 
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Figure 8. Stanton Number Distributions for 20° Sphere-Cone 
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AX [SYMMETRIC 


R n *4.0 ft (1.22m) 


1800 R (1000 K) 


V TO * 25, 000 fps (7620 mps) 0 C • 20 deg 
Altitude * 280,000 ft (85340 m) 


Re /ft • 808 
00 


Re Im * 2651 
00 


£ » 0. 1161 


Altitude - 310, 000 ft (94490 m) 

Re /ft *140.1 Re /m -460 £-0.2746 


\ 


Kang Results; 

11 Sp, TVSL, ECW 


310 Kft 


Present Results; ^ 
^ 7 Sp, TVSL, ECW 


280 Kft 


310 Kft - 


280 Kft 


Normalized Heat-Transfer Distributions for 20° Sphere-Cone 
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Figure 10. Temperature Profiles for 20° Sphere-Cone at 280 Kft 
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Figure 11. Temperature Profiles for 20° Sphere-Cone at 310 Kft 
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Figure 12. Electron Concentration Profiles for 20° Sphere-Cone at 280 Kft 


87 



, number of electrons/cm 









°(K) x 10" 


AX I SYMMETRIC 


Figure 


Altitude - 230, 000 ft (70100 m) 

Re* /ft -8630 

R n »0.5 ft (0.15 m) 

V*, » 25, 000 fps (7620 mps) 
CURVE GAS BODY THEORY 


0 C - 9 deg 

Re /m« 28314 
00 

T w - 1800 °R (1000 °K) 
ITER SHOCK WALL 


Present Results 


1 

o-o 2 

Hyp 

TVSL 

2 

NSS 

ECW 

2 

7 Sp 

Hyp 

TVSL 

2 

NSS 

ECW 

3 

CNJ 

? 

o 

Cone 

TVSL 

2 

NSS 

ECW 

4 

7 Sp 

Cone 

TVSL 

2 

NSS 

ECW 


Kang Results 



0 20 40 60 80 100 


s * s /R 


M 

n 


14. Shock-Temperature Distributions for RAM C Conditions, 230 Kft 
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Figure 15. Shock-Layer Thickness Distributions for 9° Sphere-Cone, RAM C 
Conditions, 230 Kft 
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Figure 16. Temperature Profiles for 9° Sphere-Cone Near Probe Location, RAM C 
Conditions, 230 Kft 



91 




AX I SYMMETRIC 


Altitude » 230, 000 ft (70100 m) 0 C -9 deg 

Re^ /ft -8630 Re^/m -28314 

R n - 0.5 ft (0. 15 m) T w • 1800 °R (1000 °K) 

• 25, 000 fps (7620 mps) 


CURVE 

GAS 

BODY 

THEORY 

ITER 

SHOCK 

WALL 




Present Results 



1 

7 Sp 

Hyp 

TVSL 

2 

NSS 

ECW 

2 

CVJ 

9 

o 

Hyp 

TVSL 

2 

NSS 

ECW 

3 

0-0 2 

Cone 

TVSL 

2 

NSS 

ECW 

4 

7 Sp 

Cone 

TVSL 

2 

NSS 

ECW 




Kang Results 



5 

11 sp 

Cone 

TVSL 

2 

ss 

ECW 



T° ( K ) x 10 ” 3 

Figure 17. Temperature Profiles for RAM C Conditions 
at s/R n = 90, 230 Kft 


92 




N e xlO" , number of electrons/cm 


EXPERIMENTAL DATA 


o • Fixed -Bias Probe 
□ Swept -Bias Probe 
• ♦ Affected by Probe Heating 

J Experimental Uncertainty 

Boundary-Layer Results of Evans, Schexnayder and Huber 

Present Results; 7 Sp, TVSL, ITER -2, ECW, NSS, 



Figure 18. Present Electron Concentration Profiles Compared with 

Experimental and Boundary-Layer Theory Profiles, s/R =8.8, 
RAM C Conditions, 230 Kft " 
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Figure 19. Present Electron Concentration Profiles with Different Reaction 
Rates Compared with Experimental and Kang and Dunn Profiles, 
s/R = 8.8, RAM C Conditions, 230 Kft 
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Figure 20. Stanton Number Distributions for 9° Sphere-Cone, RAM C Conditions 
230 Kft 
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Figure 21. Present Predictions of Electron Concentration Profiles without Shock 
Slip Compared with Experimental Data and Predictions of Kang and Dunn, 
s/R =8.8, RAM C Conditions 
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Figure 22. Present Predictions of Electron Concentration Profiles with Shock 

Slip Compared with Experimental Data and Predictions of Kang and Dunn, 
s/R = 8.8, RAM C Conditions 


97 


























*, * 

s = s /R 

n 


Figure 23. Normalized Surface-Pressure Distributions for 7.5° Sphere-Cone, 
Ames Conditions 
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Figure 24. Normalized Heat-Transfer Distributions for 7.5° Sphere-Cone, 
Ames Conditions 
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Figure 25. Comparison of Predicted Heat-Transfer Distributions for Sphere-Cones 
RAM C and Ames Conditions 
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